Structure and thermodynamics of platelet dispersions 
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Various properties of fluids consisting of platelike particles differ from the corresponding ones of 
fluids consisting of spherical particles because interactions between platelets depend on their mutual 
orientations. One of the main issues in this topic is to understand how structural properties of such 
fluids depend on factors such as the shape of the platelets, the size polydispersity, the orientational 
order, and the platelet number density. A statistical mechanics approach to the problem is natural 
and in the last few years there has been a lot of work on the study of properties of platelet fluids. In 
this contribution some recent theoretical developments in the field are discussed and experimental 
investigations are described. 



INTRODUCTION 



The symmetry of the shape of platelike particles is re- 
duced compared with spherical symmetry so that inter- 
actions between them depend not only on the separation 
between their centers but also on their mutual orienta- 
tions. As a result pair correlation functions and phase 
diagrams of fluids consisting of such particles differ from 
the corresponding ones of fluids consisting of spherical 
particles. In this article various properties of platelet 
fluids are discussed. The structure of the article is as 
follows. In Section II some examples of platelet fluids 
including the experimental investigations are presented. 
The knowledge about the shape of the particles and their 
size polydispersity allows one to calculate intermolecular 
pair correlation functions and the equation of state in the 
isotropic phase within the framework of an interaction 
site integral equation theory as is discussed in Sections 
III and IV. The theoretically obtained static structure 
factor and isothermal compressibility are compared with 
experimental data. Density functional theory is used in 
Section V to study orientational order in spatially homo- 
geneous bulk phases. Moreover, some brief details about 
properties of platelet fluids in contact with substrates are 
recorded. Some aspects of spatially inhomogeneous bulk 
phases are discussed in Section VI. Finally, Section VII 
draws conclusions. 



II. SHAPE OF PARTICLES AND 
POLYDISPERSITY IN SIZE 



In the following we will present examples of systems 
that illustrate the range of interest in platelet disper- 
sions. The shape of the particles and the polydispersity 
in size can be revealed by cryogenic transmission electron 
microscopy and scattering methods. 



Cryogenic transmission electron microscopy of 
polyethylene nanoplatelets 



Polyethylene is a commodity polymer that has be- 
come ubiquitous over the past several decades because 
of its low price and good mechanical properties. Hence, 
the number of applications of the material is huge and 
many millions of tons are produced worldwide annu- 
ally. However, polyethylene has hardly played any role 
in the field of nanotechnology. This is due to the prob- 
lem that polyethylene is produced either by free radical 
polymerization under high pressure and temperature or 
with metal-organic catalysts working exclusively under 
strictly water-free conditions. Polymer nanoparticles and 
their composites with inorganic compounds, however, are 
very often produced in aqueous systems. Recently, it has 
been demonstrated that ethylene can be polymerized in 
aqueous systems in a catalytic fashion by Ni(II) com- 
plexes. By virtue of this novel synthesis, long chains of 
polyethylene can be generated in a well-controlled en- 
vironment and at ambient temperature. Thus, it could 
be shown that aqueous polyethylene dispersions can be 
produced. This polymerization has opened the way of 
well-defined polyethylene nanoplatelets consisting of the 
smallest single crystal of polyethylene ever reported 
Semi-crystalline polyethylene in extremely mild condi- 
tions of pressure and temperature have been prepared. 
During this process, nanoplatelets exclusively made of 
polyethylene are formed by the chains as soon as they 
polymerize. In principle, the crystallization could start 
while the chains are still growing. However, the fast poly- 
merization process and the slow nucleation in confined 
nanoparticles prevent the crystallization from starting 
before the whole particle is formed. It is thus a melt 
crystallization that takes place in each nanoplatelet with 
an extraordinary degree of supercooling which could not 
be achieved by any other method Q. 

Specimens for cryogenic transmission electron mi- 
croscopy (cryo-TEM) have been prepared by vitrification 
of a thin liquid film of the diluted dispersion supported 
by a copper grid in liquid ethane. Figure [T] displays a 
micrograph of the particles in dilute aqueous solution as 
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obtained by cryo-TEM. The dispersion consists of flat 
platelets. The different gray scales for different particles 
can be easily rationalized by different viewing angles: If 
the normal of a platelet is nearly perpendicular to the 
electron beam, the length of the optical path through 
the particles is much longer than for parallel arrange- 
ment. The platelets have a hexagonal shape with rather 
straight edges. From cryo-TEM images a lateral dimen- 
sion (pseudo-diameter) of 25.4 (±4.3) nm and a corre- 
sponding thickness of the platelets of 6.3 (± 0.8) nm have 
been determined The observed nanoplatelets can 

be considered as single lamellae of polyethylene. The 
crystallinity of the nanoplatelets has been confirmed by 
wide-angle x-ray scattering measurements performed on 
the dispersion. The crystalline lamellae visible in Fig. Q] 
are expected to be covered by an amorphous layer as 
polyethylene is a semicrystalline polymer. However, it is 
not possible to detect this by cryo-TEM. This is due to 
the fact that the electron density of amorphous polyethy- 
lene is virtually the same as that of the low density amor- 
phous ice. Hence, there is not a sufficient contrast be- 
tween a possible amorphous layer and the surrounding 
medium. However, this problem can be treated by ap- 
plying contrast variation small-angle x-ray scattering as 
an analytical tool as is outlined in the following subsec- 
tion. Figure Q] also demonstrates that the particles are 
well-dispersed in the aqueous medium; virtually no ag- 
gregates are found. This is a prerequisite for a meaningful 
investigation by scattering methods in solution. 

Cryo-TEM has also been used to analyze the shape 
and size polydispersity of larger colloidal platelets such 
as nickel hydroxide [Ni(OH)2|_platelets of 180 nm di- 
ameter and 10 nm thickness [3|, |j| as well as gibbsite 
[Al(OH)3] platelets of 200 nm diameter and 15 nm thick- 
ness 
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However, a quantitative determination of the thickness 
of these colloidal platelets in solution is not straightfor- 
ward, because in particular thin platelets can hardly ever 
been viewed edgewise on TEM images, while very thick 
platelets touching the TEM grid with their rim have been 
found. The influence of confining surfaces on the orien- 
tation of the platelets was also discovered in suspensions 
of Ag, Ni, and Fe^Ni^ nanoplatelets [U [H, HI, [I| . In 
order to promote the edgewise orientation of platelets on 
the TEM grid one may use deliberately flocculated sam- 
ples [10( . The fact that the domains of liquid-crystalline 
phases of colloidal platelets at higher particle number 
densities extend over rather large distances is also due to 
the influence of the rigid walls of the sample cells. 



B. Scattering methods 

The theory of static scattering has been presented in 
various books. Here, we review the results necessary for 
this study. Light scattering, small-angle neutron scat- 
tering, and small-angle x-ray scattering determine the 
scattering intensity /(q, p) as a function of the scattering 
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FIG. 1: Cryogenic transmission electron microscopy micro- 
graph of polyethylene nanoplatelets in aqueous dispersion 
1]. The grey background is the low-density amorphous ice 
in which the particles are dispersed. The particles are flat 
platelets, appearing as rods when their normal is perpendic- 
ular to the electron beam (upper inset) and hexagons when 
their normal is rather parallel to the electron beam (lower 
inset). 



vector q and the number density of the dissolved particles 
p. The absolute value of the scattering vector is given by 
q = |q| = (47rn/A) sin(6'/2) in which n is the refractive 
index of the medium, A is the incident wavelength, and 
9 is the scattering angle. The scattering intensity can be 
written as 

I{^P) = Plo(<l)S(q.,p) + phncoh , (1) 

where io(q) is the scattering intensity of a single parti- 
cle, ^(q, p) is the structure factor related to the mutual 
interactions of the solute particles, and Ii nco h is the in- 
coherent contribution in small-angle neutron scattering 
that is mainly due to the hydrogen atoms of the parti- 
cles. The incoherent contribution Iincoh must be sub- 
tracted carefully from experimental data in order to ob- 
tain meaningful results on the structure and interaction 
of the dissolved particles. There is no such contribution 
in light scattering and small-angle x-ray scattering. The 
scattering intensity of a single particle may be used to 
define the form factor P(q) that is normalized to unity 
at q — according to 

J (q) = V p 2 (ApfP( q ), (2) 

where V p is the volume of a dissolved particle and Ap 
is the contrast of the solute resulting from the difference 
of the average scattering length density of the dissolved 
particles and the scattering density of the solvent. In 
the case of light scattering and small-angle x-ray scat- 
tering, the contrast is fixed by the solvent, the partial 
specific volume of the solute, and its chemical composi- 
tion. In the case of small-angle neutron scattering, Ap is 
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an available model parameter that can be varied by mix- 
ing protonated and deuterated solvent. The form factor 
of randomly oriented monodisperse circular platelets of 
radius R and thickness L is given by 



P{q) = 4 / dx 



Jl(qR\/l-x 2 ) sin 2 (qLx/2) 
(qRVT^) 2 (qLx/2) 2 



(3) 



where J\{x) denotes the cylindrical Bessel function of 
first-order. The effect of size polydispersity of the 
platelets is taken into account by an appropriate aver- 
age according to 

(7 (q)) - J dR J dLI (q)G(R,L), (4) 



where G(R, L) is a distribution function characterizing 
the degree of polydispersity. 



1. Form factor of stilbenoid dendrimers 




FIG. 2: Chemical structure of a stilbenoid dendrimer of the 
third generation (R = OCsHiz). 



Dendrimers are synthetic macromolecules with defined 
architecture that are synthesized by iterative controlled 
reaction steps. Starting from a focal unit, subsequent 
generations are connected to this initial core which re- 
sults in a treelike structure. Dendrimers composed of 
flexible units adopt a so-called dense core structure char- 
acterized by an average segmental density that has its 
maximum at the center of the molecule. This is eas- 
ily derived from the fact that flexible dendrimers can 
assume a great number of conformations in which the 
terminal groups can fold back. Hence these flexible den- 
drimers do not exhibit a well-defined surface given by 
the terminal groups of the last generation. The aver- 
age density profile thus derived can be used to calculate 
the interaction of flexible dendrimers in solution. These 
structures are well understood by now [25[. Much less 
attention has been paid to dendrimers consisting of rigid 
units [H [H [H, EH H HH, H HI, H . Figure [U shows 
the chemical structure of a rigid stilbenoid dendrimer 
of the third generation which is solely composed of stiff 
units. Starting from a central phenyl group, all subse- 
quent generations are built up by trans-stilben units, only 
the terminal groups are substituted by hexyloxy groups 
in order to ensure better solubility in common solvents 
such as toluene. Full conjugation in trans-stilbene can be 
achieved in the completely planar conformation. How- 
ever, the potential energy for a slight torsion around the 
single bonds is low in the ground state of trans-stilbene. 
Molecular modeling without taking into account the sol- 
vent has demonstrated that a stilbenoid dendrimer of the 
third generation exhibits a relatively compact platelike 
structure [28j . 

The spatial structure of the stilbenoid dendrimers in 
dilute solution has been investigated by small-angle x-ray 
scattering and small-angle neutron scattering [35j ]. Both 



methods give different information: Small-angle neutron 
scattering data taken at high contrast between the solute 
and the solvent can be used to explore the shape of the 
entire molecule. Small-angle x-ray scattering, however, is 
mainly sensitive to the stilbenoid scaffold of the molecule 
and hardly detects the hydrocarbon chains attached to 
its periphery due to the fact that these hexyl chains are 
virtually matched by the solvent toluene. 

The form factor P(q) of stilbenoid dendrimers of the 
third generation as obtained by small-angle neutron scat- 
tering is shown in Fig. [3J The form factor agrees with 
the one calculated numerically for a monodisperse cir- 
cular platelet with a radius R = 2.4 nm and thickness 
L = 1.8 nm according to Eq. ([3]). The good agreement 
between the experimental and theoretical results demon- 
strates that a full understanding of the spatial structure 
of stilbenoid dendrimers of the third generation has been 
achieved |35| . Dendrimers consisting of rigid units ex- 
hibit a rather well defined structure in solution and may 
serve as model systems for interacting monodisperse par- 
ticles in statistical physics (3(| as is discussed in Section 

El 



2. Form factor of laponite platelets 

Laponite consists of platelike clay particles with a 
radius R i=s 11.5 nm and thickness L as 0.9 nm. Each 
platelet carries a few hundred elementary charges. For 
a long time, aqueous suspensions of laponite have been 
extensively inve stig ated as a model system for platclikc 
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FIG. 3: The form factor P(q) of monodisperse stilbenoid den- 
drimers of the third generation [351 ] as obtained by small- 
angle neutron scattering (circles) and calculated according to 
Eq. ©. 



FIG. 4: The form factor P(q) of polydisperse laponite 
platelets [85| as obtained by light scattering (triangles) and 
small-angle x-ray scattering (squares). The solid line repre- 
sents the theoretical results for polydisperse circular platelets 
according to Eqs. (J2J) - ((4j) . 



theoretical description of such a well-characterized model 
clay appeared feasible. Early attempts focused on the 
structure of the electric double layer around a single 
platelet [4(| US [44[ , or around two parallel platelets [50] ■ 
Much of the recent theoretical work was initiated by con- 
sidering laponite platelets and their associated electric 
double layers as nonintersecting platelets carrying a con- 
stant electrostatic quadrupole moment [46j]. The spatial 
arrangement of these objects has been investigated by 
Monte Carlo simulations. The results of these studies 
point to a reversible sol to gel transition at platelet num- 
ber density p > 0.06 R~ 3 , where the structure of the gel 
phase is strongly reminiscent of a " house-of-cards" . 

A screened electrostatic potential between two arbi- 
trary oriented, charged platelets has been worked out 
within the framework of a linearized Poisson-Boltzmann 
theory [45|, [53|, [U, [64| . At a fixed center-to-center dis- 
tance, this repulsive electrostatic potential is maximized 
for coplanar platelets, which corresponds to the maxi- 
mum overlap of electric double layers, and minimized 
when the platelets are coaxial and parallel. An inter- 
mediate potential is found for T-shaped perpendicular 
platelets. At present only preliminary results of Monte 
Carlo simulations for platelets interacting through such 
an anisotropic potential are available, indicating that the 
platelet fluid is in the isotropic phase for platelet num- 
ber density p < 0.1 i?~ 3 [6l|. Later an interaction site 
model of laponite was put forward, and investigated by 
molecular-dynamics simulations [54j and integral equa- 
tions 57]. In this model each platelet carries discrete 
charged sites rigidly arranged on an array inscribed in 
the platelet. Sites on different platelets interact via a re- 
pulsive screened Coulomb interaction. Upon varying the 
platelet number density and screening length, a sol, gel, 



and crystal phase have been identified. Additional rim 
charges of opposite sign have been found to lead to T- 
shaped pair configurations and clustering of the platelets. 

The theoretical studies on gelation and structure for- 
mation in laponite suspensions have assumed a system 
of monodisperse platelets, i.e., thin platelets with a ra- 
dius R and a thickness L. However, it has become 
possible to visualize directly the laponite platelets by 
various techniques such as atomic force microscopy or 
electron microscopy. From these micrographs it is ob- 
vious that laponite exhibits a considerable size polydis- 
persity [HI, [H, US H3]- The degree of polydispersity 
has been quantified in a comprehensive study of laponite 
suspensions in dilute solution using light scattering and 
small-angle x-ray scattering [85]. Figure 2] displays the 
measured and calculated form factor P(q) of laponite 
platelets. From a comparison of Fig.[3]and Fig.[3]it is ob- 
vious that size polydispersity must be included in the case 
of laponite platelets because the form factor of monodis- 
perse platelets given by Eq. ((3]) has a wavy shape which 
is smoothed out by size polydispersity in the case of the 
laponite platelets. The solid line in Fig. 0] shows the re- 
sulting fit obtained for a Schulz-Zimm-distribution with 
a polydispersity expressed through R w /R n — 1.6. The 
weight-average radius R w is given by 11.5 nm, the thick- 
ness is found to be L=0.9 nm, and the number-average 
radius is denoted by R n . It is worthwhile to reiterate 
that laponite platelets exhibit an appreciable polydisper- 
sity that can be deduced from this analysis. Obviously, 
this must be taken into account when modeling the in- 
teraction of the particles at finite volume fraction. 
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III. INTERMOLECULAR PAIR 
CORRELATIONS IN THE ISOTROPIC PHASE 

Using the form factor as input into generalized 
Ornstein-Zernike equations of an interaction site integral 
equation theory, one can calculate structural properties 
of interacting platelets. In this section we will focus on re- 
sults obtained from intermolecular pair correlation func- 
tions in the isotropic phase. Calculated structure factors 
will be compared with experimental data. 



A. RISM and PRISM theory 

The systems under investigation are dispersions or so- 
lutions, but in view of the mesoscopic scale of the par- 
ticles, the solvent will be considered as a structureless 
continuum. Spatial pair correlations of an isotropic fluid 
of identical particles, each carrying n distinct interaction 
sites, are characterized by a set of intermolecular site-site 
total correlation functions hij (r, p) , where the indices i 
and j run over sites on each of two particles, and p is the 
particle number density. These functions are related to a 
set of intermolecular site-site direct correlation functions 
Cij (r, p) by the generalized Ornstein-Zernike relations of 
the "reference interaction site model" (RISM), which in 
Fourier space read [H, Wf\ 

h tj(<l,p) = ^2 UJ vm{q)c mo {q,p){uJo ] {q) + ph oj (q,p)) , 

m,o— 1 

(5) 

where the 0Jij{q) are the Fourier transforms of the in- 
tramolecular correlation functions for rigid particles. For 
flexible particles the intramolecular correlation functions 
depend on the particle number density and follow from 
a statistical average over particle configurations. In the 
case of rigid particles such as the platelets considered 
here, the 0Jij{q) are independent of the particle number 
density because the particles are not deformed due to in- 
termolecular interactions for typical concentrations in the 
fluid state. The set of generalized Ornstein-Zernike equa- 
tions must be supplemented by a set of closure relations. 
If the interaction sites are simply the centers of exclu- 
sion spheres of diameter d, to account for steric effects, 
a convenient closure is the Percus-Yevick approximation 

11,113 

hij(r, p) = — 1, r < d, Cij(r, p) = , r > d. (6) 

The experimentally accessible structure factor S(q,p) is 
defined as 



is the particle-averaged total correlation function. The 
particle-averaged intramolecular correlation function 



S(q,p) = 1 + P 



h(q,p) 
P(q) 



where 



1 " 

h(q,p) = -j T] h mo (q,p) 



(7) 



(8) 



(9) 



m,o— 1 



characterizes the geometry of the distribution of the sites, 
and hence the geometric shape of the particles. While 
the particle-averaged intramolecular correlation function 
accounts for the interference of radiation scattered from 
different parts of the same particle in a scattering ex- 
periment, the local order in the fluid is characterized by 
h(q,p) or S(q,p). 

The RISM has been proved to be a successful theory of 
the pair structure of many molecular fluids (for a review 
see Ref. [8!|). In the case of macromolecular and colloidal 
systems, with very large numbers of interaction sites, the 
number of coupled RISM equations becomes intractable, 
and a considerable simplification follows from the as- 
sumption that the direct correlation functions c%j (?, p) 
are independent of the indices i and j. This leads to 
the "polymer reference interaction site model" (PRISM) 
theory first ap plie d by Schweizer and Curro to long flex- 
ible polymers |90j. PRISM neglects end effects in that 
case. The resulting single generalized Ornstein-Zernike 
equation of the PRISM reads 

h(q, p) = P(q)c(q, p)P(q) + pc(q, p)h(q, p)P{q) , (10) 

where c(q, p) = ^ 0=1 c m „(ij, p). Taking into account 
size polydispersity leads to coupled generalized Ornstein- 
Zernike equations of the multicomponent system consist- 
ing of N species: 

h aj (q,{p a }) = P a (q)c al (q,{p a })P 1 {q) 
N 

+ P a(q) ^2 C ^(q, {Pa})psh8 7 (q, {p a }) ■ 

(11) 
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This set of N(N + l)/2 independent Ornstein-Zernike 
equations must be supplemented by as many closure re- 
lations between each pair of total and direct correlation 
functions. The PRISM integral equation theory has been 
successfully app lied to various systems, such as rodlike 
viruses jU, [92[ [9^, [94[ , polymers [H, Q , mixtures of 
spherical colloids and semiflexible polymers [97|, bottle- 
brush polymers (98|, and polyelectrolyte brushes [99] ]. 
Moreover, it has been demonstrated that the simpler 
PRISM theory yields results in good agreement with the 
more elaborate RISM calculations for lamellar colloids 

fiool. 
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B. Structure factor of stilbenoid dendrimers 

In Fig. [5] the experimental structure factors S(q,p) for 
monodisperse stilbenoid dendrimers of the third gener- 
ation (see Fig. \2§ are compared to the results of the 
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FIG. 5: Experimentally determined structure factors S(q, p) 
for monodisperse stilbenoid dendrimers of the third gener- 
ation (see Fig. [2} together with the results of the theoret- 
ical predictions of the PRISM integral equation theory ac- 
cording to Eqs. ©, 0, and JTOJl 36]. The volume frac- 
tion of the dendrimers increases from top to bottom: <j> — 
0.009,0.019,0.039,0.061,0.1. For reasons of clarity, the up- 
per data sets and lines have been shifted up by 0.1, 0.2, 0.3, 
0.4, respectively. 



integral equation theory for the PRISM [3a]. The cal- 
culated form factor P(q) (see the solid line in Fig. [3|) 
has been used as input into the generalized Ornstein- 
Zernike equation [Eq. (|10p ]. The volume fraction is given 
by 4> = V p p, where V p is the volume of an individual parti- 
cle. The generalized Ornstein-Zcrnikc equation has been 
solved numerically together with the Percus-Yevick clo- 
sure. From Fig. [5] it is apparent that the PRISM integral 
equation theory is rather accurate. The magnitude and 
the scattering vector range of the suppression of S(q, p), 
i.e., the deviations from the value 1 at small scattering 
vectors, are characteristic for the size and the shape of 
the dendrimers as well as the volume fraction. On the 
basis of the experience with both PRISM and RISM it is 
expected that the results of the integral equation theory 
for the RISM would lead to very similar results, provided 
the same form factor is used [100| . 



C. Scattering intensity of polyethylene 
nanoparticles 

The above mentioned dispersions consisting of polydis- 
perse polyethylene nanoplatelets (see Fig. [I} contained 
in addition just enough of the surfactant sodium dodecyl 
sulfate to stabilize the particles against coagulation 
The surface tension of 65 mN/m of the aqueous disper- 
sion demonstrated that virtually all the surfactant was 
adsorbed onto the particles leading to an effective charge 
of the platelets which has to be taken into account in 



the PRISM integral equation theory. One may consider 
a multicomponent system involving N species of charged 
platelets with number densities p a , where 1 < a < N. 
Each platelet of species a contains n a equivalent inter- 
action sites. The interaction potential between sites on 
particles of species a and 7, carrying the charges z a e and 
z a e, is of the generic form: 



U ay (r) 



er 



■ exp(— Kpr) 



(12) 



where A^> = is the usual Debye screening length. 
The system is an aqueous dispersion, but in view of the 
mesoscopic scale of the platelets, the solvent is mod- 
elled as a structureless dielectric continuum providing a 
macroscopic permittivity e. Any microscopic counteri- 
ons or ions from added electrolyte will be considered at 
the linear response (or Debye-Hiickel) level, i.e., they will 
screen the electrostatic potential due to the interaction 
sites on the charged platelets on a scale given by the De- 
bye screening length. The underlying assumption entails 
that the charge distribution on the mesoscopic particles 
does not contribute to screening. 

The Laria-Wu-Chandler closure relation [1011 ] has been 
used in order to supplement the Ornstein-Zernike equa- 
tions [Eq. (jlip ] of the multicomponent system: 

h ay (r, {p a }) = ln[h aj (r, {p a }) + 1] 

+ P a * [c Q7 + (k B T)~ 1 n a n J U al ] * P 7 (r) , 

(13) 

where the asterisk * denotes a convolution product. In 
addition, steric effects have been taken into account using 
the Percus-Yevick approximation. The cryo-TEM micro- 
graph displayed in Fig. [1] shows that the particles are 
hexagonal platelets with a pseudo-radius of 12.5 nm and 
a thickness of 6.3 nm. As argued above, however, the 
cryo-TEM micrograph shows only the crystalline lamella. 
Since polyethylene is a semicrystalline polymer, an amor- 
phous layer in which the chains fold back must exist on 
both sides of the lamella. Hence, there must be a thin 
amorphous layer on both sides of the particles that needs 
to be taken into account. The form factor of such a "ham- 
burger" is given by 



dx 



JiiqRVT 



(qRVT 



r 2\2 



L 8iniqLx/2) Ap, 



sin( g W2) 
(qL c x/2) V lc n} 

x (LAp! + L c (Ap c - Api)y 2 



(qLx/2) 
2 



(14) 



Here R denotes the radius of the platelet, L is the overall 
thickness of the particles, while L c denotes the thickness 
of the crystalline lamella. Api and Ap c arc the contrasts 
of the outer two platelets and of the inner platelet, respec- 
tively. The polydispersity of the particles estimated from 
the cryo-TEM micrographs has been taken into account 
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in the fitting using a Gaussian distribution in radius and 
thickness. The closed set of equations (|TTj) - (fT4"|) has been 
solved numerically by a standard iterative procedure in 
order to obtain the scattering intensity |l[ . The solid lines 
in Fig. [6] display the form factors. The small-angle x-ray 
scattering intensities of the polyethylene nanoplatelets 
have been measured at six different contrasts starting 
from a stock solution of (/> = 0.017 of the nanoplatelets 
dispersed in pure water. The different contrasts have 
been adjusted by adding different amounts of sucrose. 
The volume fractions of added sucrose increased from 
0.0 (lowest contrast) to 0.378 (highest contrast) while 
the corresponding volume fraction of the nanoparticles 
decreased from tfi = 0.017 to cf> = 0.007. For the sake of 
clarity the scattering intensities related to different con- 
trast have been shifted vertically in Fig. The figure 
demonstrates that varying the contrast leads to marked 
differences in the scattering intensities. In particular, 
the maxima of the scattering intensities are shifted in a 
characteristic manner when changing the contrast. The 
overall dimensions following from the theoretical descrip- 
tion are the weight-average radius R w = 12 nm and the 
weight-average thickness L w — 9 nm. The thickness of 
the crystalline layer is L c = 6.3 nm. From the particle 
volume, the polyethylene density and the chain molecular 
weight, one can estimate that each particle is made up of 
ca. 14 chains. Moreover, from the entire series of small- 
angle x-ray scattering-intensities the contrast results to 
Api = 309 electrons/nm 3 for the amorphous part and 
Api = 333 electrons/nm 3 for the crystalline part. The 
electron densities of the crystalline and amorphous parts 
of the particles can be used to understand the cryo-TEM 
micrograph in detail. The value for Api is close to the 
electron density of low density amorphous-ice. Therefore 
one can conclude that the amorphous parts are virtually 
matched in the cryo-TEM. 

For small magnitudes of the scattering vector q the cal- 
culated scattering intensities for noninteracting particles 
(solid lines) on the one hand, and the integral equation 
results for interacting particles (dotted lines) as well as 
the experimental data (symbols) on the other hand de- 
viate due to strong repulsive interactions between the 
particles. The observed decrease of the scattering inten- 
sity at small scattering vectors is considerably more pro- 
nounced than for fluids consisting of hard platelets such 
as the rigid dendrimers discussed above. From this result 
one can conclude that the colloidal stability is achieved 
by electrostatic repulsion between the platelets brought 
about by the adsorbed surfactant. It has been estimated 
that each polyethylene nanoplatelet is surrounded by a 
layer consisting of a few hundred surfactant molecules [l| . 



D. Structure factor of laponite 

While the preceding subsections have demonstrated 
that the PRISM theory can describe pair correlations 
of both uncharged and uniformly charged platelike parti- 
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FIG. 6: Measured scattering intensity of polyethylene 
nanoplatelets as a function of the magnitude of the scattering 
vector q (symbols) 0]. All intensities have been normalized to 
the volume fraction <j> of the platelets in the dispersion. The 
volume fraction of the added sucrose increases from bottom 
to top (0, 6.2, 10.3, 18.0, 25.4, 37.8 vol.%) while the volume 
fraction cj> of the nanoplatelets decreases from bottom to top 
(1.73, 1.56, 1.44, 1.23, 1.02, 0.68 vol.%). The five lowermost 
intensities are shifted down by a factor of 10, 10 2 , 10 3 , 10 4 , 
10 , respectively. The solid lines represent the result of the 
modeling of the small-angle x-ray scattering data assuming a 
dispersion of noninteracting polydisperse platelets. The short 
dashed lines (q < 0.25 nm -1 ) represent the scattering in- 
tensity calculated for a dispersion of interacting platelets as 
obtained from the PRISM integral equation theory [Eqs. (|11[) 
- I|14p]. The differences between the dotted and solid lines re- 
flect the intermolecular interaction between the nanoplatelets. 



cles in solution, disagreement between experimental and 
theoretical results has been found in the case of laponite 
platelets if a purely repulsive screened Coulomb poten- 
tial is used in t he t heoretical calculations [Hj]. It has 
been suggested [l02j | that at low concentrations the ef- 
fective interaction between laponite platelets is charac- 
terized by a competition of long-range Coulomb repul- 
sion and short-range van der Waals attr action sim i lar t o 
colloidal systems and protein solutions jl03l . Il04 1 1051 ] . 
Moreover, the charge density is higher on the face than 
on the rim of a platelet, which leads to a modification of 
the interaction potential between sites on particles within 
the computationally demanding multicomponent PRISM 
theory. Therefore, a mesoscopic coarse graining, whereby 
particles act via an effective potential U(r) has been pro- 
posed |85| . To this end the Ornstein-Zernike equation 

%> p) = c{q, p) + pc(q)h(q, p) (15) 

has been solved together with the hypernetted-chain clo- 
sure relation 

h(r, p) = In[ft(r, p) + l] + c{r, p) + (ksT^Uir) . (16) 
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FIG. 7: The structure factor S(q, p) of laponite platelets for 
various volume fractions (open circles: 0.02 %; open squares: 
0.04 %; open triangles: 0.08 %; solid circles: 0.12 %; solid 
squares: 0.16 %). The solid lines represent the results of the 
integral equation theory [Eqs. (|15|) - (|17p ] using the effective 
potential Z7(r) shown in the inset [85l ]. 



These equations follow from Eqs. (fTTj) and (fT5)) for AT = 1 
and P a (q) — 1. The calculated structure factor 



<S(g,p) = 1 + ph(q,p) 



(17) 



is compared in Fig. [7] with experimental data [85]. For 
all volume fractions, good agreement with the experimen- 
tal results has been achieved with a density-independent 
effective potential U(r) that is attractive for small dis- 
tances but repulsive for larger distances (see the inset in 
Fig. [7J . The effective potential fulfills the sta bility con- 
dition J °° drr 2 U(r) > discussed by Ruelle [l06| | and 
provides a representation of the underlying many-body 
interactions in the system. In general, to any S(q, p) 
there corresponds a unique effective center of mass po- 
tential U e ff(r,p), capable of reproducing S(q,p), irre- 
spective of the underlying many-body interactions in the 
system. This effective potential is rather independent of 
the particle number density for the laponite suspensions, 
i.e.,U e f f (r,p) «U(r). 

Since structure factors and pair correlation func- 
tions of polyelectrolyte solutions calculated within the 
PRISM framework and using a purely repulsive screened 
Coulomb potential are in good agreement with experi- 
mental data and computer simulations, it is worthwhile 
to compare scattering data for polyelectrolytes with the 
data obtained for laponite. Figure [8] displays S(q = 0,p) 
together with S(q = q max , p) for various polyelectrolytes 
and laponite [851 ] . where q max is the absolute value of 
the scattering vector at the main peak. The struc- 
ture factors of the polyelectrolytes exhibit a main peak 
S(q = qmaxiP) > 1 anci a small value S(q — 0, p) < 0.11 
for small scattering vectors, while the structure factors of 
laponite are characterized by a main peak and a rather 
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FIG. 8: The structure factor at zero scattering vector 
S(q — 0, p) (lower symbols) together with the structure factor 
at the main peak S(q = q m ax,p) (upper symbols) for vari- 
ous suspensions: laponite at five volume fractions: 0.02 %, 
0.04 %, 0.08 %, 0.12 %, 0.16 % from left to right (squares 
[8^]): polystyrenesulfonate of length L = 6.8 nm at two con- 
centratio ns: .1 mol/L and 0.2 mol/L from left to right (up 
triangles [l07j ]): polystyrenesulfonate of length L = 40 nm 
at two concentrations: 0.1 mol/L and 0.2 mol/L from left to 
right (down triangles |l07ll): DNA of length L — 57 nm at 
0.05 mol/L (circles [l07j l: proteoglycan at two different salt 
concentr ation s: 0.05 mM and mM from left to right (plus- 
symbols EM): DNA of length L = 380 nm at 0.05 mol/L 
(diamonds [107[ |): tobacco mosaic virus at five concentrations: 
0.11 mg/ml, 0.27 mg/mL, 0.43 mg/mL, 1.05 mg/ml, 2.07 
mg/mL from left to right (crosses |l09l ]). 



large value 0.7 < S(q = 0, p) < 0.95 for small scattering 
vectors. The enhancement of the small angle (q — > 0) 
value of S(q, p) of laponite as compared to the polyelec- 
trolytes signals increased density fluctuations. Since the 
polyelectrolyte solutions remain liquid like even at higher 
concentrations opposite to the laponite suspensions, the 
observed qualitative different behavior of the structure 
factors of laponite may be considered as indicative of the 
sol to gel transition at higher concentrations. In addi- 
tion the values of S(q = 0, p) for the liquid like laponite 
suspensions are smaller than 1 and decrease with increas- 
ing concentration in contrast to strong small q upturns 
which have been observed experimentally for clay suspen- 
sions in the gel phase, polyelectrolyte gels and mixtures, 
and low ionic strength polyelectrolyte solutions. These 
strong upturns signal strong concentration fluctuations 
indicative of aggregation, or spinodal instability reminis- 
cent of the behavior observed in computer simulations 
and PRSIM integral equation theory of mixtures of op- 
positely charged particles [97j ■ 
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IV. EQUATION OF STATE AND ISOTHERMAL 
COMPRESSIBILITY IN THE ISOTROPIC PHASE 

The structure factor provides a direct link with ther- 
modynamics via the compressibility equation [Hf 



1 + 2B 2 p + 3B 3 p 2 + 0(p 3 



(21) 



lim S(q, p) = pk B TK T {p) , 



(18) 



where kt(p) is the isothermal compressibility. The os- 
motic pressure P(p) (equation of state) then follows from 



P(p) 
k R T 



dp'S- 1 (q = 0,p'). 



(19) 



Various attempts have been made to develop accurate 
theories for the equation of state of fluids consisting of 
nonspherical particles: 

(a) Scaled particle theory, which is very successful for 
hard sphere fluids, has been exte nded to prolate and 
oblate ellipsoids of revolution |lldj |. however with mod- 
erate succ ess when gauged against Monte Car lo simula- 
tions 111]. Recently, it has be en shown ll 12il that the 
results of scaled particle theory |ll2l . Ill3l Ill4j for plate- 
like particles or the closely related model of hard cut 
spheres are in disagre e ment with co mputer simulation 
data 0, [HI, EH El Ell ITT1IT20] . 

(b) Onsager theory jl2ll Il22j |. based on the second 
virial coefficient alone, can be "rescaled" Il23l . Il24l Il25| . 



Although this semi-empirical procedure leads to reason- 
ably good results for rodlike par ticle s, it is much less 
satisfactory for platelike particles [l26j | . 

(c) Many theoretical studies on hard sphere fluids and 
depl etion agents use the so-called free volume theory 
[l27| | , in which the free volume accessible to a single par- 
ticle plays a major role. This free volume theor y ha s 
been studied within a fundamental measure theory [Till ]. 
However, it has been demonstrated that the resulting 
third virial coefficients of the equation of state for both 
hard cylinders and hard cut spheres differ from com- 
puter simulation results (see tables I and II in Ref. [112| |). 
Theoretical approaches based on fundamental measure 
theory do not yield correct third virial coefficients and 
equation of states due to the occurrence of so-called 
lost cases, i.e., the fact that configurations of three par- 
ticles with pairwise overlap but no triple overlap do 
not c ontribute to t hermodynamic properties (see, e.g., 

Refs. [m EH EH). 

These theoretical and computer simulation studies 
demonstrate that the understanding of thermodynamic 
properties of nonspherical particles needs to be improved. 
Recently, the measured inverse structure factor S^ 1 (q = 
0, p) extrapolated to vanishing scattering vectors of stil- 
benoid dendrimers of the third generation h as been m od- 
elled in terms of the so-called y3 expansion jl3ll . Il32j 



S- 1 (q = 0,p) = 
1 + 2{B 2 V P - 2)cf> - 



(3B 3 V*-8B 2 V p + 6)<t> 2 



(1 



(20) 



The ?/3 theory reproduces the exact second and third 
virial coefficients, B 2 and -B3, respectively. However, its 
practical applicability is limited due to the difficult nu- 
merical evaluation of the third vi rial c oefficient B3 in 
the case of nonspherical particles jl33l |. while the sec- 
ond virial coefficient B 2 for an iso tropi c hard convex 
body fluid is known exactly (see Ref. 134] and references 
therein): 



B 2 = V„ 



ApRp . 



(22) 



Here A p and R p — (l/47r) J dA p H p are the surface 
area and the mean radius, respectively, where the lo- 
cal mean curvature is denoted as H p . For a circular 
platelet of radius R and thickness L one has V p — ttR 2 L, 
A p = 2ttR(R + L), and R p = irR/4 + L/4. For platelike 
stilbenoid dendrimers of the third generation the second 
virial coefficient B 2 /V p — 5.54 as calculated from Eq. (|22|) 
with R = 2.4 nm and L = 1.8 nm agrees with the exper- 
imentally determined second virial coefficient [36| . 

In the framework of both scaled particle theory and 
Rosenfeld's fundamental measure theory S^ 1 (q = 0,p) 
is also given by Eqs. (|20[) and (f2Tj) but the third virial 
coefficient reads 



B 



(SPT) 



V p 2 + 2R P A P V P + ^R 2 p A 2 p (23) 



within scaled particle theory and 



B 



(FMT) 



v 2 



2R p A p Vp 



1 A3 

12tA 



(24) 



within fundamental measure theory. In Fig. [5] 
the experimentally determined inverse structure factor 
S^ 1 (q — 0,p) extrapolated to vanishing scattering vec- 
tors of stilbenoid dendrimers of the third generation is 
compared with the results of scaled particle theory ac- 
cording to Eqs. (HOI), and (HJ) and fundamental 
measure theory according to Eqs. (f20|) . (|22|) . and (|24)l . 
With increasing volume fraction the theoretical results 
of both scaled particle theory (dashed line) and funda- 
mental measure theory (dotted line) deviate from the ex- 
perimental data (symbols). These deviations are mainly 
due to the fact that the predicted third virial coefficients 
B ( 3 SPT) and B ( 3 FMT) are too small. However, the PRISM 
theory discussed above leads to a agreement with the ex- 
perimental data as is apparent from Fig. [5] [36] . 

This comparison of experimentally determined 
S^ 1 (q — 0, p) of platelike particles with the predictions 
of the well-known scaled particle and fundamental 
measure theory has confirmed earlier caveats con- 
cerning the applicability of th ese t heories to freely 
rotating nonspherical particles [l 12| . For example, 
b {fmt) = 2^2^6/3 = 6 580i? 6 ag ca i cu i atec i f rom 

Eq. (|2"4")l for thin circular platelets in the limit L — > 
underestimates the third virial coefficient Ufa = 10.8 29-R 6 
according to computer simulation studies [lid . IllE 
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FIG. 10: Illustrations of the polar angle 9 and the azimuthal 
angle of the normal a; of a platelet of radius R and thickness 
L. The center of mass of the particle is located at r = (0, 0, 0). 



0.1 



0.2 



FIG. 9: Inverse structure factor S~ 1 (q = 0,p) extrapolated 
to vanishing scattering vectors of stilbenoid dendrimers of the 
third generation (with the same symbol code as in Fig. 
The dashed line follows from the scaled particle theory ac- 
cording to Eqs. I|20")l. and while the dotted line 
represents the results of the fundamental measure theory as 
obtained from Eqs. (T22]), and (JUJ) [3f|. 



V. ORIENTATIONAL ORDER IN SPATIALLY 
HOMOGENEOUS BULK PHASES 

In this section we will outline some points on the basic 
description of orientational order in platelet dispersions 
and we will record details about the phase behavior of 
spatially homogeneous bulk fluids. 



A. Density functional theory 



intrinsic Helmholtz free energy of an ideal gas of nonin- 
teracting particles, 

N 

F id [{ Pi }} = fc B T^ / d 3 rdu Pi (r,u>) 

x (ln[47rA? Pi (r,w)]-l) , (26) 

where Aj is the corresponding thermal de Broglie wave- 
length. The excess (over the ideal gas) free energy func- 
tional -Fe£c[{Pi}] i s m general a very complicated, highly 
non-trivial object, because it is a characteristic prop- 
erty of an interacting many-body system. The functional 
^e£c[{pi}] i s dealt with in various ways, which specify the 
explicit forms of the theory. 

By means of a diagrammatic expansion it can be shown 
that the exact excess free energy functional in the low 
density limit reduces to a virial approximation, which is 
the starting point of the following considerations. Within 
a third-order virial approximation the excess free energy 
functional is given by 



The systems under investigation are dispersions involv- 
ing N species of colloidal particles. The number density 
of the centers of mass of platelets of species i at a point 
r with an orientation u> = (9, <j>) of the normal of the 
platelets (see Fig. [TU]) is denoted by pi(r, u)) while for 
spherical structureless particles the center of mass num- 
ber density is independent of oj. The equilibrium density 
profiles of the fluid under the influence of external po- 
tentials Uj(r, minimize the grand potential functional 
plll35l . [l36t 

ti[{pi}] = F ld [{ Pl }} + F ex [{ Pl }} 

N 

+ 22 d 3 rdu}pi(r,u)(Ui(r,<jj)- Hi) .{2b) 

z=l 

Here pi is the chemical potential of particles of species i 
and Fid[{ P i}} + F ex [{ P i}] is the intrinsic Helmholtz free 
energy of the fluid (i.e., the Helmholtz free energy in 
the absence of the external potentials). -FidtlPi}] is the 



Fex[{Pi}] = 




x/ij(ri,wi;r2,w 2 )Pj(r2,W2) 



x(l + 3 d3r 3^3/j7c(r2,W2;r 3 ,u)2) 
Xp fe (r3,w 3 )/ fci (r3,a73;ri,a>i)) , (27) 

where /y(ri,u>i; r2, 0J2) is the Mayer function of the pair 
interaction potential between two particles of species i 
and j. For sterically stabilized colloidal particles the 
most important feature of their pair interaction is the 
strong mutual repulsion whenever the distance between 
the centers of mass of the particles is smaller than that at 
which their surfaces touch each other. For nonspherical 
particles this shape-dependent repulsion is a function of 
the mutual orientation of the particles. In view of the 
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steepness of the repulsive potential between two parti- 
cles, the latter may be regarded as hard particles, such 
that the pair interaction potential is infinite if the par- 
ticle volumes overlap and is zero otherwise. By match- 
ing the indices of refraction of the colloidal particles and 
the solvent it is possible to effectively switch off the om- 
nipresent dispersion forces and to create colloidal suspen- 
sions in which the actual effective interaction between the 
colloidal particles very closely resembles this hard core 
potential. In this way theoretical models of hard body 
fluids have their actual experimental counterparts with 
which their properties can be compared quantitatively. 
For such hard body fluids the corresponding Mayer func- 
tion equals —1 if the hard particles overlap and is zero 
otherwise [88]. Two thin platelets (R ^> L, see Fig. ITO]) . 
separated by a distance ri 2 , intersect if the inequality 



r 12 • (u>i X w 2 )| < jR 2 sin j 12 - (r 12 ■ u>i) 2 



i? 2 sin 2 7l2 -(r 12 -a; 2 ) 2 (28) 



with ojj = (sin 6j cos <j>j , sin 9 j sin (f>j , cos 6j ) , j = 1,2, 
and T12 = i"i — r 2 is fulfilled jll6l |. 712 is the angle be- 
tween the normals u)\ and ui 2 of the two platelets. For 
thin platelets the intersection volume is like a line seg- 
ment for uix 7^ us 2 . Configurations with exactly equal 
orientations u)% — u> 2 are not considered for freely ro- 
tating thin platelets in three dimensions because they 
have a vanishing statistical weight. Minimization of O 
with respect to pi(r,ui) leads to integral equations for 
the equilibrium density profiles: 



P»(r,w) 



1 



4ttA? 



exp 



1 dF[{ Pi }] 
k B Td Pl {r,u>) 



(29) 



These equations can be solved numerically for given 
chemical potentials pi and given external potentials 
Ui(r,uj) using a Picard scheme with retardation, i.e., 
the equations pi(r, u>) = J[{pi}} are solved according 
to Pi {to, (jj) = mix J[{pf + (1 — vnix)p { \ (r, u)) 
with a mixing parameter < mix < 1. Starting from 
an initial guess for p' ' (r, u>), the solutions of the equa- 
tions are obtained by an iterative approach I = 0,1, 2, ... . 
The mixing parameter is adjusted empirically to en- 
sure convergence. In practice, first the Mayer func- 
tions are calculated and are stored for all required val- 
ues of (ri,u>i;r 2 ,a; 2 ). Thereafter the integral equations 
[Eq. ([29])] are solved. 

For a homogeneous and isotropic bulk fluid with 
Ui(v,u)) = in a macroscopic volume V the grand po- 
tential functional [Eq. (|2"5"]) ] reduces to 



N 



Vk B T 



5> (\n(Af Pl )-l -iM(kBT)- 1 ) 

N ( . N \ 

i,j=l \ k=l ) 



(30) 



where pi = V 1 J d 3 r du: pi(r,ui) are the total particle 

number densities. B^ is the second virial coefficient 
between particles of species i and j: 



1 



2(4tt) 



d 3 rx duii du) 2 /y(ri, Uf, 0, w 2 ) . 



_ (31) 

In accordance with the standard definition L 88j , the third 
virial coefficient involving three particles of species i, j, 
k is given by 



B 



(ijk) 



1 



d 7"i d 3 r 2 du)\ du> 2 du)j, 



3(4vr) 3 

x/i J -(ri,a;i;0 1 a>2)/jfc(ri - r 2 , w 2 ; 0, u> 3 ) 



x/fei(r 2 ,a;3;0,a;i) 



(32) 



The equation of state derived from the grand potential 
[Eq. ([30]) ] takes the following form: 



1 fen 



P 



k B T k B T\dV 



t,{M 



(33) 



N N N 

?( 2 )„.„.^ \' R(3), 



\ i— 1 — 1 ij,k—l 

(34) 

The same equation without the second and the third term 
on the right-hand side holds for the ideal gas limit, i.e., 
nonintcracting particles. 



B. Fundamental measure theory for thin hard 
platelets in three dimensions 

Within a geometry-based density functional theory 
the excess free energy functional for monodisperse thin 
platelets in three dimensions is obtained by integrating 
over an excess free energy density, 

F ex [p] = £pS / d 3 r dwi du> 2 $({n„}) , (35) 

lD7T^ J 

where the spatial and angular arguments of the weighted 
densities n v are suppressed in th e notation. The free 
energy density <!> is given by |l37l ] 

$({»v}) = ni(r,u)i)n 2 (r,wi) 

+ 7^— n 2 (r, u>i)n 3 (r, a>i; u> 2 )n 2 (r, w 2 ).(36) 

Z47T 

The weighted densities are related to the number density 
of the center of mass of the platelets p(r, u)) according to 

/duj* 
—wx{r,u)'\u} X )*p{T,u}'), (37) 

«2(r,u>i)= W2(r,u)i)*p(r,u>i), (38) 

/du)' 
—w 3 (r,uj ;w 1 ;w 2 )*p(r,w ),(39) 
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where the asterisk ★ denotes the spatial convolution and 
the weight functions are given by 

i«o(r,«i) =i(5(fl-|r|)<y(r-a;i) ) (40) 



«i 2 (r,wi) =2Q{R-\t\)6(T-u Jl ) [ 



(41) 



u;i(r,o>i;w 2 ) =- s |a>i • (w 2 X r)|too(r,a>i), (42) 
K 
g 

U)3(r,a>i;w 2 ;<»;3) = — |a?i • (u> 2 x a; 3 )|ii; 2 (r, a>i).(43) 



Here 0(x) is the Heaviside step function and S(x) is the 
Dirac distribution. The Mayer function of the pair inter- 
action potential between two hard platelets [Eq. (j2"5|) ] is 
recovered via the convolution 

/(r,u;i;0,u> 2 ) = -2w 2 (r, wi) * wi(r, on; u> 2 ) . (44) 

However, the fundamental measure theory does not yield 
the correct third virial coefficient due to the above men- 
tioned occurrence of lost cases [137j . The fundamental 
measure approximation for the third virial coefficient is 
nonvanishing only for three platelets at position vectors 
ri,r 2 ,r3 with a common triple intersection. The lost 
cases are compensated by overcounting the cases with 
triple intersections, such that the resulting equation of 
state is in reasonable agreement with the exact results 
for the isotropic bulk fluid. In the case of thin platelets 
this has been achieved by using the prefactor l/(247r) in 
the second term in Eq. (|36[) which leads to the third virial 



coefficient B 



(FMT) 



2n 2 R 6 /3. The equation of state as 



obtained from this fundamental measure theory is indeed 
in reasonable agreement with computer simulation data 
as is apparent from Fig. 1111 

The density functional theory allows one to consider 
the possibility of an isotropic to nematic phase transi- 
tion in a monodisperse fluid of platelets as a function of 
particle number density. The broken symmetry in the 
spatially homogeneous nematic phase may be character- 
ized by the nematic order parameter 



S=— !— [ d 3 rdujp( r,uj)( 3 cos 2 9 - 1) 
8npV J 



(45) 



which varies between in the isotropic phase and 1 when 
all platelets are oriented along one axis. For the bulk 
isotropic (I) to nematic (N) phase transition the follow- 
ing values of the particle number density and the nematic 
order parameter at two-phase coexistence have been ob- 
tained from fundamental measure theory: p/i? 3 = 0.42, 
p N R 3 = 0.46, 5/ = 0.05, and S N = 0.49 [l37|. These 
values are similar to the ones obtained from computer 
simulations while a second virial approximation based 
on Eq. (f2"7| without the last term overestimates both 
the density jump and the nematic order parameter at 
coexistence and yields piR 3 — 0.67, pnR 3 = 0.85, and 
Sn = 0.79. Even for thin platelets the chance that one 
platelet has contact with two or more platelets simultane- 
ously is significant, and thus for this fluid the third virial 




FIG. 11: Equation of state of thin platelets (R S> L, see 
Fig. I10p as obtained from fundamen tal m easure theory accord- 
ing to Eqs. (|35fl - (1431) ( solid line [l37j l ) a nd from computer 
simulations (circles [1161 ] and squares [46j]). The horizontal 
line is the tie-line illustrating isotropic-nematic two-phase co- 
existence according to the fundamental measure theory. 



coefficient is non- negligible. For comparison we note that 
the contribution of the third virial coefficient to the equa- 
tion of state of a fluid consisting of thin rods is vanishing 
small, reflecting the smal l pro bability that three rods will 
simultaneously intersect [122j |. 

Density functional theory has been applied to fluids 
consisting of th in hard platelets in contact with a hard 
substrate [138]. Platelets lying very close to the sub- 
strate must adopt nearly a fully parallel alignment due 
to interactions with the substrate. The probability of 
finding platelets touching the substrate is increased com- 
pared with the bulk due to effective entropic interactions. 
Upon increasing the bulk platelet den sity t he substrate is 
completely wetted by a nematic film [l39j . For the fluid 
confined by two parallel hard substrates, a first-order 
capillary nematization transition for large slit widths has 
been determined [l40j |. 

The enrichment and ordering of clay platelets near 
surfa ces is important for the oil and gas production 
14l|| . The first stage of recovering oil and gas from 



an underground reservoir is the drilling of wells into 
the hydrocarbon-bearing rock. This is performed by 
rotary drilling using, for example, water-based drilling 
fluids that contain smectic clay platelets such as ben- 
tonite or montmorillonite. Among the many functions 
that drilling fluids must possess to be successful for a 
drilling operation are: carrying the drilling cuttings and 
transporting them back to the surface; suspending the 
drilled cuttings when circulation is stopped; cooling and 
cleaning the bit; reducing the friction between the drilling 
string and the sides of the hole; preventing the inflow 
of fluids from the permeable rocks that are penetrated; 
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maintaining the stability of the wellbore; and forming a 
thin, low permeability filter-cake which seals the pores 
and other openings in formations penetrated by the bit. 
Scanning electron microscope images of mixed-metal hy- 
droxide - bentonite filter-cake formed on the surface of a 
sandstone rock show in detail how the mixed-metal hy- 
droxide and bentonite platelets form bridges across the 
entrance of rock pores of size 30 /jm and larger, although 
the r adius of an individual platelet is less than 0.5 /im 
[l4l| . Effective entropic interactions together with van 
der Waals interactions at submicron distances play an 
important role in understanding the underlying mecha- 
nisms of orientational ordering near surfaces [142| . Den- 
sity functional theory provides the basic de scrip tion of 
such inhomogeneous colloidal platelet fluids [134]. 



5. 




C. Strictly two-dimensional hard platelets 

The density functional theory presented in the previ- 
ous subsection is based on the weight functions given 
by Eqs. (|40| - (|43|) and the deconvolution of the Mayer 
function [Eq. (|44|) ]. Although the particle shape is equiv- 
alent, overlapping pair configurations are very different 
for freely rotating platelets in three dimension as com- 
pared to the strictly two-dimensional case. While for thin 
platelets in three dimensions the intersection volume is 
like a line segment (ui =^ W2), the intersection volume 
is an area for platelets in two dimensions (ui = u^). 
Rosenfeld's weight funct ions for two-dimen sional hard 
platelets of radius R are pM EM EM EH 

w 2 (r) = 0(iJ-|r|), (46) 
«;i(r) = S(R-\r\), (47) 

wi(r) = «U2(r)|^, (48) 
w (r) = Wl (r)^. (49) 

The exact Mayer function is approximated by 

/(r) = 2w (r)*w 2 (r) 

+ — (wi (r) * wi (r) - wi (r) * wi (r)) , (50) 
Ztt 

where r is the center to center vector and the asterisk * 
denotes the two-dimensional convo lutio n. The evaluation 
of the integrals in Eq. (JSDJ) yields jl37l | 

/» = 

e(2i? " r) (^ arccos (^) + wiF = ^)- (51) 

Figure [T2] illustrates that /(r) and /(r) differ qualita- 
tively. In particular, /(r) exhibits a divergence in the 
case that the two platelets touch each other, i.e., for |r| = 
2R. However the deviations from the exact value (solid 



r/(2R) 

FIG. 12: The exact Mayer function f(r) for platelets of radius 
R in two dimensions (solid line) and the approximate repre- 
sentation f(r) according to the fundamental measure theory 
(dashed line) as obtained from Eq. (|5ip . respectively. The 
distance between the centers of the platelets is denoted by r. 
Note that f(r) diverges upon approaching r = 2R. 



line in Fig. [T2")) balance such that the second virial coef- 
ficient is exact, i.e., J d 2 r /(r) = J d 2 r /(r). Higher or- 
der virial coefficients have been calcu lated recently using 
molecular dynamics simulations |l46| . Moreover, the fun- 
damental measure theory for two-dimensional platelets 
has been considered as a useful tool to study inhomoge- 
neous systems such as laser in duced freezing and melting 
of confined colloidal particles |l47l ] despite the deficien- 
cies of the theory in reproducing the Mayer function. 



D. Colloidal mixtures of hard spheres and platelets 

Depletion interactions between large colloidal particles 
induced by smaller particles, which can be either solvent 
particles or a collodial component in its own right, are of 
significant research interest because of the importance of 
these effective interactions in various colloidal processes. 
For example, flocculation of colloids can be driven by 
the addition of non-adsorbing polymers via the depletion 
mechanism. On the other hand, stabilization of colloidal 
suspensions can be achieved by the addition of adsorbing 
polymers. Whereas experimental and theoretical studies 
have focused on binary hard-sphere fluids as well as on 
colloidal mixtures of hard spheres and hard rods or poly- 
mers, less atte ntion has been paid to hard p l atelets ac ting 
as depletants [111, EM EM EM EM EM EM E3 • 

Density functional theory for thin hard platelets has 
been used to investigate the depletion potential between 
two hard spheres due to the presence of thin platelets 
[l50j |. The calculated depletion potential between two 
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spheres exhibits an attractive primary minimum at con- 
tact when the spheres touch each other. This mini- 
mum deepens upon increasing the platelet density and 
a small repulsive barrier at larger sphere separations 
develops. However, the depletion barrier is typically 
less than the thermal energy fcgT, and therefore un- 
likely to significantly alter the kinetics of aggregation of 
the hard spheres at platelet densities smaller than one 
half the density of the isotropic phase at bulk isotropic- 
nematic two-phase coexistence. Nonetheless, with in- 
creasing platelet density the integrated strength of the 
effective interaction between the spheres becomes signifi- 
cantly weaker and thus reduces the thermodynamic onset 
of flocculation [l50j |. 

A fundamental measure-based density functional the- 
ory has been used to study binary colloidal fluids con- 
sisting of hard spheres and thin platelets in their bulk 
and near a planar hard substrate. Figure [T3] displays 
the calculated phase diagram for a binary mixture of 
spheres and thin platelets for a size rati o Rs/R = 2, 
where R s is the radius of the spheres jl53l | . The tie-lines 
are horizontal because the coexisting phases are charac- 
terized by the equality of the chemical potential of the 
platelets /i in the coexisting phases. The binodal for 
coexisting states is shown, for which a sphere-rich and 
a platelet-poor liquid phase coexists with a sphere-poor 
and a platelet-rich liquid phase. The coexistence region 
is bounded by a lower critical point below which only 
a single stable phase is found. Upon decreasing \x the 
adsorption behavior of the mixture near a planar sub- 
strate changes qualitatively so that two cases have to be 
distinguished. For > —5.551 the sphere-rich 

phase does not wet the wall at two-phase coexistence. 
The layer thickness of the sphere-rich phase forming close 
to the wall increases continuously upon approaching the 
bulk phase boundary, but remains finite at coexistence. 
For ^/(fcfiT) < —5.551, however, the wall is wetted com- 
pletely. The transition to complete wetting is first-order 
because the excess adsorptions of both the spheres and 
platelets along the bulk coexistence curve jump to macro- 
scopic values upon crossing the wetting transition point. 

Binary mixtures of silica spheres of radius 
R s = 350 nm and silica-coated gibbsite platelets of 
radius R = 229 nm and tickness L = 56 nm have been 
inves tigated by means of confocal microscopy and TEM 
|l54l | . The depletion- induced aggregation of the spherical 
particles has been observed in a TEM micrograph. In 
contrast to sphere-polymer and sphere-rod mixtures, 
the addition of platelets to a solution of spheres close to 
the fluid-solid binodal of the latter slowed down early 
crystalline ordering of the spheres, which led to more 
grain boundaries and powder-like structures. It has been 
argued that in the two-phase region the disordered sedi- 
ment structures are kinetically arrested, which prevents 
further crystallization. Moreover, fluid-like microphases 
of platelets were trapped in the sediment which is due to 
the simultaneous sedimentation of platelets and spheres 
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FIG. 13: Bulk and surface phase diagrams of binary mixtures 
of spheres of radius R 3 and thin platelets of radius R = R s /2 
as a function of the chemical potential o f the platelets fi and 
the number density of the spheres p s jl53l |. The straight 
dashed lines are tie-lines illustrating liquid-liquid phase co- 
existence. The solid and open circles denote the bulk critical 
point and the wetting transition point, respectively. Between 
the wetting transition point and the critical point the sphere- 
rich liquid phase completely wets the interface between a hard 
substrate and the sphere-poor liquid phase. 



as well as due to depletion interactions. 

It has been shown that spheres and platelets can be 
sepa rated by adding a third species to the suspension 
[l55| . The aggregation of platelets (R — 0.5 ^m and 
L = 0.2 ^m) and spheres (R s = 0.5 /im) of 1-eicosene 
wax in a dilute aqueous suspension has been induced 
by nanometer-sized spherical micelles of sodium dode- 
cyl sulfate. Below the critical micelle concentration nei- 
ther the platelets nor the spheres aggregate. Upon in- 
creasing the micelle concentration the platelets aggregate 
into columns due to effective entropic attraction while 
the spheres remain unaggregated. The columnar aggre- 
gates of platelets resemble the stacked-coin structures of 
red blood cells, known as rouleaux, which form when 
polymers are added to blood. Analogous to fractioned 
crystallization, the aggregated platelets can be repeat- 
edly separated from the unaggregated spheres by gravita- 
tional creaming, yielding a shape-specific colloidal purifi- 
cation method. For high micelle concentrations not only 
the platelets aggregate but also the spheres aggregate 
with other spheres. It has been argued that the observed 
shape-selective colloidal separation might be applicable 
to a broad class of mixtures of spherical and nonspheri- 
cal c olloids such as viruses, organelles, proteins, and cells 
PI - 

Very recently it has been investigated how the surface 
roughness of colloidal pl atelet s can affect the depletion 
attraction between them [l56l |. The depletion attraction 
between polymeric pentagonal platelets (pseudo-radius 
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FIG. 14: Within the Zwanzig model, platelets are represented 
by square parallelepipeds of thickness L and width D, which 
can take only three orientations, along the x, y or z directions. 
The number density of the center of mass of the platelet in 
the figure is denoted by p x {r). 



R w 1.5 pirn and thickness L = 1 /zm) can be suppressed 
when the nanoscale surface asperity heights on the face 
of the platelets become larger than the size of nanometer- 
sized spherical micelles of sodium dodecyl sulfate acting 
as depletion agent. In the opposite limit, the attraction 
reappears and columnar stacks of platelets are formed. 
Exploiting this, the site-specific roughness on only one 
side of the platelets has been increased in order to form 
aligned platelet dimers. The two smoother surfaces of 
two Janus platelets aggregate face-to-face, exposing the 
rougher surfaces that cannot aggregate. The general 
problem of how surface roughness affects the strength 
and range of the depletion attraction is interesting and 
needs to be elucidated in future. 



E. Density functional theory for the Zwanzig model 

Taking into account the simultaneous presence of ori- 
entational ordering and polydispersity in density func- 
tional calculations for freely rotating nonspherical parti- 
cles is computationally demanding because the integrals 
to be evaluated in Eq. ([77]) are high-dimensional, already 
within a third-order virial approximation. In order to re- 
duce this computational effort one ma y use the Zwanzig 
model for hard square parallelepipeds jl57l |. Within the 
Zwanzig model particles of species i = 1, ...,N are repre- 
sented by rectangular blocks of size Li x D, x Z)j. The 
positions of the centers of mass vary continuously, while 
the orientations of the normal of each particle are re- 
stricted to discrete directions f3 — x,y,z (see Fig. IT4"|) . 

Using the notation ot£ (x, z) — oti(x, z,8 = 7r/2, <f> = 0), 

oty (x,z) = oti{x,z,6 = tt/2,4 1 = 7r /2), and ai l \x,z) = 
ai{x,z,6 = 0, 4> = 0) for a = p,U, the grand potential 
functional [Eq. ([25)) ] can be written as 

n[{pfm = F ld [{pf{v)}]+F ex [{pf{v)}] 



+ E E j <?rpf(v)(uf(v)-p) 

i=l fi=x,y,z 

(52) 

with the ideal gas contribution [Eq. ([26]) ] 

Fu[{ P f{v)}] = 

N 

fc ^ T EE / d3 rpf(r) (MAfpjftr)] - l) .(53) 

The third-order virial approximation of the excess free 
energy functional [Eq. (|27p] reads 

F ex [{pf(r)}} = 

N 

-^f E E fM^^f^^p^) 

x(l + lEE/^/S(r 2 ,r3) 

xp^)/^ (i^rO), (54) 

where /^ J X(ri,r2) is the Mayer function which for hard 
bodies as considered here equals —1 if two particles of 
species i and j with orientations f3\ and fj% overlap and 
is zero otherwise. With the definition 

s % = A + {Li - D t )S a(3 , (55) 

which represents the spatial extent in direction a = 
x, y, z of a particle of species i with orientation (3 of the 
normal, the Mayer function can be written explicitly as 

-U @ (l{ S «k+ S «l)-\r»i-r a 2\) , (56) 

where r a \ is the projection of the position vector ri 
in a direction. The density functional theory is com- 
pletely specified by the excess free energy functional and 
the Mayer functions. Minimization of f2 with respect 

to p^(r) leads to integral equations for the equilibrium 
density profiles. 

Zwanzig's model may be considered as a coarse-grained 
version of the Onsager mod el, which allows for contin- 
uously varying orientations [l2ll Il22j |. The model of- 
fers the advantage that the determination of phase dia- 
grams and density profiles becomes numerically straight- 
forward, allowing one to study various aspects of platelet 
fluids in detail and to scan wide parameter ranges. More- 
over, Rosenfeld's very successful fundamental measure 
theory for hard sphere systems has been extended to 
the Zwanzig model, which is an important alternative in 
view of the uncertain convergence behavior of the virial 
serie s of t he equation of state of platelets of finite thick- 
ness [l58j . Within the fundamental measure theory, one 
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postulates the follo wing form of the excess free energy 
functional [HI EM [l6l[ : 

F ex [{pf(r)}} = k B T J d 3 r d>({p«(r)}) , (57) 
with the reduced excess free energy density 

$ = -n ln(l - n 3 ) + h — , (58) 



1 - n 3 (1 - n 3 y 



and the weighted densities 

JV 



n <( r ) = EE / ^^(r-nJpW^) (59) 

i=l /3 ^ 

for Z e {0, lx, ly, lz, 2x, 2y, 2z, 3}. The weight functions 
u>i l p (r) are obtained by expressing the Fourier transform 
of the Mayer function [Eq. ([56]) ] as a sum of products of 
single particle functions. The explicit expressions of the 
weight functions read: 



<4« = I II s ( s «l/ 2 -\ r e\) ■ (60) 

a— x.y,z 

= II e(^S/2-M) , (61) 



w 26 5«/2-|r Q | ( ) 

W la/?( r ) = 



*(^/2-|r a |) 



2e(S»/2-|r Q | 



<M r )> ( 62 ) 
W $(r). (63) 



In the limit of monodisperse thin platelets (L <C D, see 
Fig. [T4]) the excess free energy functional [Eqs. (|57|) - 
()63|) ] reduces to a third-order virial approximation. It is 
worthwhile to emphasize that the fundamental measure 
theory for the Zwanzig model yields the exact second 
and third virial coefficients irrespective of the shape and 
size polydispersity of the particles. For monodisperse 
platelets the second virial coefficient is given by Eq. 
with V p = LD 2 , A p = 2D(D + L), and R p = (L + 2D)/3, 
while the third virial coefficient reads 



B (ZWA) 



V p 2 + 2R P A P V P + 2A% . 



(64) 



As a direct application of the fundamental measure 
theory one may consider the possibility of an isotopic 
to nematic phase transition in a monodisperse fluid of 
platelets, as a function of the aspect ratio L/D. The 
broken symmetry in the nematic phase may be charac- 
terized by the nematic order parameter [Eq. (|45p ]: 



1 



2pV 



d 3 r {2 Pz (v) - p x (v) - Py {r)) , (65) 



where p = V 1 J d 3 r (p x (r) + p y (r) + p z (r)). In the ne- 
matic phase S = 1 and p x = p y . The one-component 
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FIG. 15: Equation of state of a monodisperse fluid consisting 
of rectangular Zwanzig platelets of surface size D x D and 
thickn ess L (see Fig.fT?]and Eqs. (I52|l - (|63p ) for various aspect 
ratios [158]: L/D = 0.01,0.1,0.2,0.3 (from bottom to top). 
The horizontal lines are tie-lines illustrating isotopic-nematic 
two-phase coexistence. 



version of the grand potential [Eq. ((52)) ] has been min- 
imized wi th re spect to the order parameter S for each 
density p 158]. The resulting equations of state in the 
isotropic (S — 0) and nematic (S ^ 0) phases are plotted 
versus p in Fig. 1151 together with the horizontal tie-lines 
connecting the two phases, for several aspect ratios L/D. 
The value of the effective packing fractions 4>i = piD 3 
and 4>n = PnD 3 of the coexisting phases are seen to be 
rather insensitive to the aspect ratio L/D, while the first- 
order phase transition narrows as L/D increases, i.e., 
A0 = <pN — <j>j decreases with increasing L/D. The re- 
duced pressure PD 3 /(ksT), however, increases rapidly 
with L/D, as one might expect. 

One may next consider the richer case of highly asym- 
metric binary mixtures of large and small platelets. Poly- 
dispersity, both in diameters and in thickness, is intrinsic 
to the aforementioned experimental gibbsite samples. In 
the case of mixtures of spherical colloidal particles large 
size asymmetry may lead to depletion-induced phase sep- 
aration. In the case of platelets depletio n- ind uced segre- 
gation competes with nematic ordering [162], as is also 
the case for mixtures of long and short rods, or of thin 
and thick rods. Figure [TBI displays a phase diagram of bi- 
nary mixtures of square platelets constructed as function 
of the chemical potential u 2 of the small platelets and the 
number density p\ of the large platelets [D\ > D%). The 
tie-lines are horizontal because of equality of u 2 of the co- 
existing phases. For large negative chemical potential p 2 , 
the systems exhibit a first-oder isotropic (I) to nematic 
(N) phase transition and the density gap at the isotropic 
to nematic phase transition broadens with increasing u 2 
in agreement with calculations for binary mixtures of in- 
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FIG. 16: Phase diagram of a fluid consisting of a binary 
mixture of large thin Zwanzig platelets [L\/D\ = 0.01] and 
small Zwanzig platelets [D1/D2 = 2, L2/D2 = 0.16] accord- 
ing to Eqs. (|52|l - (|63|) as a function of the chemical poten- 
tial of the small platelets (12 and the density of the large 
platelets pi [l58t ] ■ The straight solid lines are tie- lines illus- 
trating isotropic-nematic (I-N) and nematic-nematic (N1-N2) 
two-phase coexistence, respectively. The chemical potential 
at I-N coexistence of a monodisperse fluid consisting of the 
small platelets is ^/(ksT) = 0.68. 



finitely thin platelets [1631 ]. At high densities in the ne- 
matic phase, the small platelets are nearly excluded and 
hence the large platelets are in equilibrium with a reser- 
voir of small platelets. This implies that the pressure 
of the nematic phase is equal to that of a reservoir of 
small platelets at a chemical potential ^2 ■ The isotropic- 
nematic two-phase coexistence regime is bounded by an 
upper critical point above which a single stable nematic 
phase is found. The nematic phase demixes into two ne- 
matic phases (Ni, N2) at sufficiently large values of p,2- 

Upon increasing the aspect ratio L2/D2 of the small 
platelets at fixed L\ and D\, the lower critical point of 
the nematic-nematic two-phase coexistence region shifts 
to smaller values of p.2, until the nematic-nematic and 
isotropic-nematic two-phase coexistence regions start to 
overlap, giving rise to a triple point at which two nem atic 
phases (Ni, N2) coexist with an isotropic phase [l58j | . 

Figure [T7l displays an alternative representation of the 
phase diagram shown in Fig. \W\ in terms of the num- 
ber densities of the large and the small platelets p\ and 
P2, respectively. The figure illustrates how the composi- 
tion of the mixture varies upon increasing the chemical 
potential (and hence the reservoir density) of the small 
platelets. Whereas the density of the thin platelets p\ 
is always smaller in the isotropic phase than in the ne- 
matic phase along the coexistence curve, a density inver- 
sion for the thick platelets has been found. The density 
of the thick platelets P2 is larger in the isotropic phase 
than in the coexisting nematic phase for small chemi- 



FIG. 17: Phase diagram of a binary platelet fluid, correspond- 
ing to Fig. [16] in the density-density, i.e., pi-p% plane [l58| ]. 
The dotted lines indicate phase boundaries. The chemical 
potential of the small thick platelets is kept fixed for each 
solid line and increases from bottom to top: p,z/(kBT) = 
-5,-4,-3,-2,-1,0,1,1.9,3. 



cal potentials /i2, while P2 is smaller in the isotropic 
phase than in the coexisting nematic phase for large 
chemical potentials p.2- This remarkable phenomenon 
of isotropic-nematic density inversion has been observed 
experimentally for__a_ dispersion of sterically stabilized 
gibbsite platelets 



111 ]. The calculated density profile 



P2{z) of thick platelets at t he isotropic-nematic inter- 
face exhibits oscillations jl64| . The interfacial tension for 
platelets with their mean bulk orientation parallel to the 
interface is smaller than for the corresponding perpen- 
dicular configuration. The calculated interfacial tension 
7 « 10 nN/m fo r gib bsite platelets agrees with the mea- 
sured value [l(| Il65| . Both the large platelet and small 
platelet density profiles exhibit pronounced oscillations 
on one side of the nematic-nematic interface, provided 
the chemical potentials are sufficiently high. Upon ap- 
proaching the lower critical point of the nematic-nematic 
demixing region (see Fig. [TO]) the oscillations vanish at a 
Fisher- Widom line and the interface broadens 



A second-order virial approximation [Eq. (|54[) . without 
the last term] has been used to investigate the phase dia - 
gram of binary mixtures of rods and platelets [l66l . Il67l ] . 
It has been found that mixing rods and platelets can sta- 
bilize a biaxial nematic phase in which the symmetry axes 
of particles of different species point along mutually per- 
pendicular directions. This phase behavior was later con- 
firmed in calculations for a polydisperse mixture of r ods 
and platelets using fundamental measure theory jl68l ] . A 
virial approximation has also been used to study wet- 
ting and capillary nematization of platelet and rod fluids 
in cont act with p lanar and geometrically structured sub- 
strates jl63l Il69l ] . The predicted complete wetting of the 
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substrate-isotropic liquid interface by a uniaxial nematic 
platelet film has been found later in a colloidal disper- 
sion of sterically stabilized gibbsite platelets 165]. The 
wetting behavior of the nematic gibbsite phase has been 
revealed by polarized light microscopy. 

The relaxation dynamics of Zwanzig platelets has been 
inve stiga ted by means of dynamic density functional the- 
ory [170| | . The driving force for time dependence has been 
expressed in terms of gradients of the local chemical po- 
tential. The formalism has been applied to model an ini- 
tially homogeneous stable or metastable isotropic fluid 
which is perturbed by switching a two-dimensional array 
of Gaussian laser beams. Switching on the laser beams 
leads to an accumulation of platelets in the beam centers. 
If the initial chemical potential and the laser power are 
large enough, a preferred orientation of the particles oc- 
curs, breaking the symmetry of the laser potential. After 
switching off the laser beams again, the system can follow 
different relaxation paths. It either relaxes back to the 
homogeneous isotropic state or it forms an approximately 
elliptical high-density core which is elongated perpendic- 
ular to the dominating orientation in order to minimize 
the surface free energy. For large supersaturations of the 
initial isotropic fluid, the high-density cores of neighbor- 
ing laser beams of the two-dimensional array merge into 
complex superstructures. 

Recently, nonequilibrium steady states in an open sys- 
tem connecting two reservoirs of platelike colloidal par- 
ticles have been investigated using a d ynam ic density 
functional theory for the Zwanzig model [l7l| . Inhomo- 
geneities of the local chemical potential generate a dif- 
fusion current which relaxes to a nonvanishing value if 
the two reservoirs coupled to the system sustain differ- 
ent chemical potentials. The relaxation process of initial 
states towards the steady state turned out to comprise 
two regimes: a smoothening of initial steplike structures 
followed by an ultimate relaxation of the slowest diffusive 
mode. The position of a nonequilibrium interface and the 
particle current of steady states depends nontrivially on 
the structure of the reservoirs due to the coupling be- 
tween translational and orientational degrees of freedom 
of the fluid. 



F. Equilibrium properties of charged Zwanzig 
platelets 

Bulk properties and free interfaces of mixtures of 
charged platelets and salt have been studied withi n th e 
density functional theory for the Zwanzig model jl72l |. 
The excess free energy functional has been decomposed 
into the fundamental measure contribution [Eq. (|57|)] and 
an electrostatic part which has been derived by functional 
integration of an extension of the Debye-Huckel pair dis- 
tribution function with respect to the interaction poten- 
tial. The calculations have shown that the isotropic and 
nematic binodals are shifted to higher platelet densities 
upon increasing the platelet charge. The Donnan poten- 



tial between the coexisting isotropic and nematic phases 
has been inferred from bulk structure calculations. Non- 
monotonic density and nematic order parameter profiles 
have been found at a free interface interpolating between 
the coexisting isotropic and nematic phases. Moreover, 
electrically charged layers form at the free interface lead- 
ing to monotonically varying electrostatic potential pro- 
files. For fixed salt density, the interfacial tension de- 
creases upon increasing the macroion charge. 

Analytically and numerically calculated surface phase 
diagrams of charged Zwanzig platelets and salt in contact 
with a charged planar substrate exhibit first-order wet- 
ting for sufficiently small platelet charges and isotropic 
bulk order as well as first-order drying for su fficie ntly 
large platelet charges and nematic bulk order jl73| |. A 
crossover from monotonic to nonmonotonic electrostatic 
potential profiles upon varying the surface charge density 
has been observed. Nonmonotonic electrostatic potential 
profiles are equivalent to the occurrence of charge inver- 
sion. Due to the presence of both the Coulomb interac- 
tions and the hard-core repulsions, the surface potential 
and the surface charge do not vanish simultaneously, i.e., 
the point of zero charge and the isoelectric point of the 
surface do not coincide. 



VI. SPATIALLY INHOMOGENEOUS BULK 
PHASES 



The Zwanzig model offers the advantage that the 
difficult determination of spatially i nhomogeneous bulk 
phases becomes numerically feasible |l64l . Il74 Il75l Il76j . 
The predicted phase behavior is in qualitative agree- 
ment with the results obtained for free ly ro tating paral- 
lelepipeds using computer simulations 1771 ]. In order to 
find spatially inhomogeneous bulk phases such as colum- 
nar and smectic phases, the spatially homogeneous so- 
lutions of the Euler-Lagrange equations for the funda- 
mental measure theory of the Zwanzig model [Eqs. (I52|) - 
(f6"5| ] are perturbed by adding a narrow Lorentz peak and 
iterated using a Picard scheme. If the latter converges 
towards the spatially homogeneous solution (isotropic or 
nematic phase), one concludes that this is indeed the 
equilibrium state because it is stable under density vari- 
ations. This method is the numerical analog of a bifur- 
cation analysis since a narrow Lorentz peak has a broad 
Fourier spectrum. In the case of the phase diagram for 
the binary mixtures of large thin and small thick platelets 
shown in Figs. [H] and [T7| it has been found that the 
nematic phase Ni undergoes a phase transition into a 
columnar phase, where parallel columns of small platelets 
are surrounded by single platel ets w ith orientations per- 
pendicular to the column axis |l64j . For small values of 
pi, the nematic to columnar phase transition is weakly 
first-order. The nematic-columnar boundary line changes 
to strongly first-order at the tricritical point where the 
nematic- nematic boundary line ends. 

Over the past few years, spatially homogeneous and 
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inhomogeneous bulk phases of the aforementioned sus- 
pensions of colloidal gibbsite platelets have been investi- 
gated. For example, van der Kooij et al. Q have stud- 
ied two systems which differ by the degree of polydis- 
persity of the radius of the platelets (17 % and 25 %, 
respectively). The phase behavior of the platelet suspen- 
sions as function of the volume fraction (f> of the platelets 
is shown in Fig. 1181 Isotropic-nematic two-phase coex- 
istence has been observed just below cj) = 0.2, yield- 
ing an isotropic upper phase and a birefringent nematic 
bottom phase in the tubes containing the suspensions. 
Macroscopic phase separation has been found to be com- 
plete within 12 hours. The width of the biphasic gap 
A0 = <pN — 4>i increases with increasing the degree of 
polydispersity of the platelet radius. Upon increasing 
the volume fraction of the platelets to <f> « 0.4, both sus- 
pensions enter a biphasic region where a nematic phase 
coexists with a more concentrated birefringent columnar 
phase. A comparison of the experimentally observed vol- 
ume fractions at the isotropic to nematic and nematic 
to columnar phase transitions with computer simulation 
data for mondisperse platelets [1171 ] has shown that the 
transitions in the experiment are shifted to lower volume 
fractions which might be due the polydispersity in size 
and shape of the gibbsite platelets [§] . The stability of 
the observed columnar phase is remarkable in view of 
predictions for the terminal size polydispersity of about 
10 % for the crystal phase of hard spheres and the termi- 
nal length polydispersity of about 18 % for the smectic 
phase of hard rods. Small-angle x-ray scattering has re- 
vealed a he xago nal intercolumnar ordering in the colum- 
nar phase [1781 ]. The system forms a powder consisting 
of true long-range-ordered columnar crystallites, where 
the relatively large free space between the columns al- 
lows for the accommodation of rather highly polydisperse 
platelets. However, under the influence of gravitational 
compression very little space is available. The geomet- 
rical frustration induced by the platelet polydispersity 
can suppress the ordering upon increasing density and 
favor hexaticlike structuring. The resulting novel colum- 
nar phase has been found at the bottom of a tube con- 
taining a gibbsite suspension [l78j . The hexatic phase is 
characterized by short-range translational order while its 
bond-orientational order is long ranged. 



Moreover, polydisperse colloidal gibbsite platelets have 
been reported to form an opal-like columnar crystal with 
iridescent Bragg reflections [20j. It has been shown that 
the formation process of the iridescent phase under slow 
sedimentation can be accelerated by orders of magnitude 
using a centrifugation without arresting the system in 
a disordered glassy state. It has been argued that the 
formation of the hexatic phase with a finite position cor- 
relation length may play a role in the self-organization of 
the platelets [201 ]. 
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FIG. 18: Phase diagram of two gibbsite suspensions which 
differ by the degree of polydispersity of the platelet radius 
(triangles, 17 % and circles, 25 %) [Sj. The relative volume 
fraction of the nematic (N) and columnar (C) phase are plot- 
ted as function of the platelet volume fraction <f>. For low 
volume fraction the isotropic phase (I) is stable. The dotted 
lines mark the boundaries of the two-phase coexistence re- 
gions of the system with 17 % polydispersity of the platelet 
radius. 



VII. SUMMARY AND OUTLOOK 

This work has been devoted to the investigation of 
the structure and thermodynamics of platelet disper- 
sions. The shape of the particles and the polydispersity 
in size can be revealed by cryogenic transmission elec- 
tron microscopy [Fig. [Tj and the form factor [Eqs. (Q]) 

- Q] measured by light scattering, small-angle neu- 
tron scattering, or small-angle x-ray scattering [Figs. [3] 
and 0] . The combination of these methods has demon- 
strated that well-defined and perfectly dispersed single 
lamellar nanoplatelets of the cheapest polyolefin, namely 
polyethylene, can be made in water [l|. Due to their un- 
usual size, structure, and thermal history, these platelets 
will contribute to fundamental and long-standing issues 
of polymer crystallization [2j. 

The interaction site integral equation theory [Eqs. ([5]) 

- (1131) ] is quite successful in describing experimental 
structure factors of fluids containing uncharged platelets 
[Fig. \5\ or uniformly charged platelets [Fig. [B]. Rigid 
dendrimers [Fig. [2] may serve as model systems for inter- 
acting monodisperse particles in statistical physics [36(, 
while the colloidal platelets synthesized so far are highly 
polydisperse in size. 

The inverse structure factor extrapolated to vanish- 
ing scattering vectors as predicted by both scaled par- 
ticle theory and fundamental measure theory is system- 
atically smaller than the experimental data for platelike 
stilbenoid dendrimers of the third generation and the pre- 
diction of the interaction site integral equation theory 
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[Fig. [9] and Eqs. (TO]) - ([24])]. The substantial differences 
observed between the experimental data and the results 
of scaled particle theory and fundamental measure the- 
ory confirm earlier caveats concerning the applicability 
of th ese theories to freely rotating nonspherical particles 
[l 12j | and are due to the fact that both theories do not 
yield the correct third virial coefficient. Hence there is 
a clear need to improve both scaled particle theory and 
fundamental measure theory for platelets. 

The different charge densities on the rim and the face of 
laponite platelets may lead the significant difference be- 
tween the effective pair potentials and structure factors 
of laponite and polyelectrolytes [Figs. [7] and [5]. How- 
ever, at the present stage there is no definite answer to 
this issue because no theoretical calculations are avail- 
able yet. An integral equation approach based on the 
reference interaction site model with different charges on 
the face and the rim of the platelets is expected to be 
helpful for the understanding of the microscopic struc- 
ture of laponite suspensions. An important point that is 
common to many experimental studies of laponite sus- 
pensions is the analysis of obvious noncquilibrium states 
reached after passing the gel line. As a matter of fact, 
considerable efforts have been devoted to the investiga- 
tion of the exact position of the gel line. The experimen- 
tal difficulties of the investigation of such time-dependent 
states are at hand. Moreover, no clear information on 
the nature of the interparticle interaction potential can 
be derived from essentially noncquilibrium states. 

Density functional theory [Eqs. (|25|) - ([34]) ] has been 
used to study the phase behavior of thin hard platelets 
[Fig. [TO] in three dimensions. The system exhibits a first- 
order isotropic to nematic phase transition with a small 
density jump at two-phase coexistence and a rather low 
nematic order parameter in the nematic phase. Within 
a geometry-based density functional theory the excess 
free energy is expressed in terms of weighted densities 
[Eqs. (|35p - l|43p]. However, this fundamental measure 
theory does not yield the correct third virial coefficient 
due to the occurrence of lost cases. These lost cases can 
be compensated by overcounting the cases with triple 
intersection such that the resulting equation of state in 
both the isotropic and the nematic phase as well as the 
location of the isotropic to nematic phase transition agree 
with computer simulation data [Fig. [11] . 

Density functional theory has also been used to study 
binary colloidal fluids consisting of hard spheres and thin 
platelets in their bulk and near a planar hard substrate. 
This system exhibits liquid-liquid coexisting of a phase 
that is rich in spheres (poor in platelets) and a phase that 
is poor in spheres (rich in spheres) [Fig [TO] . An interest- 
ing issue that has not been addressed theoretically is how 
the r oughness of platelets modifies depletion interactions 

Although the particle shape is equivalent, overlap- 
ping pair configurations are different for freely rotat- 
ing platelets in three dimension as compared to the 
strictly two-dimensional case. Fundamental measure 



theory yields only an approximate expression for the 
Mayer function of strictly two-dimensional hard platelets 
[Fig. [TO] . The understanding of the deconvolution of the 
Mayer functi on fo r platelets in two dimensions needs to 
be improved |l79| |. 

Zwanzig's model of parallelepipeds with only three al- 
lowed orientations [Fig. [TO] has been generalized to the 
case of monodisperse, bidisperse, and polydisperse sys- 
tems of platelets. In view of the uncertainties concerning 
the convergence of the virial series [Eqs. ([52]) - {55])], a 
fundamental measure theory has been used to derive a 
free energy functional [Eqs. ([57]) - ([63]) ]. This funda- 
mental measure theory for the Zwanzig model yields the 
exact second and third virial coefficients irrespective of 
the shape and size polydispersity of the particles. In the 
monodisperse case, the isotropic to nematic phase transi- 
tion narrows with increasing aspect ratio of the platelets 
[Fig. [TO]. The bulk phase diagram for a representative 
example of binary hard platelet fluids [Figs. [TO] and [TO] 
involve an isotropic and one or two nematic phases of 
different concentrations as well as a columnar phase. 

The Zwanzig model offers the advantage that the 
difficult determination of spatially inhomogeneous bulk 
phases and density profiles becomes numerically feasible. 
The calculated low interfacial tensions are in agreement 
with experimental data for gibbsite platelets. Moreover, 
the predicted wetting behavior of platelet fluids in con- 
tact with substrates has been confirmed experimentally 
[l65j | . It is appealing to base the construction of theories 
for the dynamics of platelets on the equilibrium funda- 
mental measure functional for the Zwanzig model [l70j |. 
In addition to the technological relevance of this issue 
(e.g., sediment transport processes in hydraulic engineer- 
ing), the investigation of dynamic properties of platelet 
fluids provides the opportunity to study the coupling of 
translational and rotational motion as well as the time 
dependence of the orientation of particles induced by ex- 
ternal fields. 

The liquid-crystalline phase behavior of colloidal gibb- 
site platelets has been studied over the past few years. It 
has been observed that suspensions of both sterically and 
charge-stabilized gibbsite platelets at sufficiently high 
concentrations display iridescence. From small-angle x- 
ray scattering experiments it follows that these iridescent 
phases have a columnar structure [Fig. [TO] with hexago- 
nal intercolumnar ordering. This is quite remarkable be- 
cause the gibbsite platelets have a polydispersity of 19 % 
in diameter and 27 % in thickness. In general, size poly- 
dispersity of colloidal particles is expected to suppress 
the formation of ordered phases. 

For the case of monodisperse hard platelets an accu- 
rate description of a columnar liquid crystal phase at 
high packing fractions has been pre sented using an im- 
proved free volume theory jl80L \l8l\ . It has been shown 
that orientational entropy of the platelets in the one- 
dimensional columns leads to a different high-density 
pressure compared to the prediction from traditional cell 
theory. Quantitative agreement has been found with 
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Monte Carlo simulation results for various thermody- 
namic properties. In view of this progress made in de- 
veloping an accurate descriptions of a columnar liquid 
crystal consisting of monodisperse hard platelets, future 
work may focus on the understanding of the influence of 
platelet size polydispersity, gravity, and external fields 
on the phase behavior. Discotic liquid crystal science 
is still a rather young research field compared with the 
more developed field of liquid crystals consisting of rod- 



like molecules. Nevertheless, many applications for dis- 
cotic liquid systems in areas such as photovoltaic and 
optical compensation films for displays have been inves- 
tigated over the last years. However, the microscopic 
understanding of the properties of columnar phases and 
the alignment at surfaces needs to be improved in partic- 
ular in view of the emerging research area of smart low- 
dimensional building blocks consisting of discotic liquid 
crystals. 
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